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' This article addresses linear hyperbolic partial differential equations and pseudodifferential equa- 

tions with strongly singular coefficients and data, modelled as members of algebras of generalised 
functions. We employ the recently developed theory of generalised Fourier integral operators to con- 
struct parametrices for the solutions and to describe propagation of singularities in this setting. As 
required tools, the construction of generalised solutions to eikonal and transport equations is given 
^ ■ and results on the microlocal regularity of the kernels of generalised Fourier integral operators are 

obtained. 
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Introduction 



This article is part of a series of papers aiming at developing microlocal analysis in Colombeau algebras of 
' generalised functions [9, 10, 11, 12, 13, 20, 24, 28]. While these algebras originally have been introduced 

I for handUng nonlinear partial differential equations with distributional data, they have turned out to be 

. equally useful for treating linear hyperbolic equations with strongly singular coefficients and data in the 

past decade. Consider the Cauchy problem for the linear hyperbolic system 

m 

dtu ^^Aj{x,t)djU + B{x,t)u + F{x,t), {x,t) e M"+\ 
(0-1) j=i 

u(x,0)=5(x), xeM". 

In case the initial data g or the driving term F are distributions and the coefhcient matrices are, e.g., 
discontinuous, the system formally contains multiplicative products of distributions. Simple examples 
show [21] that such systems may fail to have solutions in the sense of distributions (referred to as classical 



* Supported by JRF, Imperial College London 

t Partially supported by FWF (Austria), grant Y237 



1 



solutions in the sequel). For second order strictly hyperbolic equations, the minimal regularity of the 
principal coefficients is log-Lipschitz regularity in time for admitting unique classical solutions [3, 4]. 
For transport cqiiations, one may go down to BV-cocfficicnts and mcasTirc data (sec the recent surveys 
[1, 16]). Differentiating the equation and trying to establish an equation for the higher derivatives of the 
formal solution already brings one beyond the admissible non-regularity, as will, e.g., letting B be a delta 
potential or _F a space-time white noise. 

The Colombeau algebra tJ(M'"+^) is a differential algebra containing the space of distributions 7?'(]R™+^) 
as a subspace and thus it provides a framework in which all operations arising in (0.1) are meaning- 
ful. Indeed, existence and uniqueness of solutions to problem (0.1) in the Colombeau algebra (as well 
as its dual) have been proven under various conditions: Aj and B real valued n x n matrices with en- 
tries in C/(M™+^), Aj symmetric [22]; symmetric hyperbolic systems of pseudodifferential operators with 
Colombeau symbols [18, 28], strictly hyperbolic systems of pseudodifferential operators [14]. 

In the classical setting, Fourier integral operators (FIOs) and microlocal methods are used to study the 
dependence of the solution on the initial data, in particular, propagation of singularities. Indeed, it 
is the purpose of this paper to derive a generalised FlO-representation of the Colombeau solution and 
to predict its (generalised) wave front set. In the Colombeau framework, pseudodifferential operators 
and their apphcation to microlocal regularity and hypoellipticity have been developed in [9, 11, 12, 
18, 20, 24], Fourier integral operators with generalised phase and amplitude have been developed in 
[10, 13]. The present paper initiates the application of generalised FIOs to hyperbolic equations with 
Colombeau coefficients. We shall show how to construct FlO-parametrices for transport equations and 
scalar hyperbolic pseudodifferential equations. Further, the generalised wave front sets of the kernels of 
generalised FIOs are determined, as well as the propagation of the generalised wave front sets through 
application of a generalised FIO. Finally, we study the generalised Hamlitonian flow and prove various 
results on propagation of singularities. 

In short, it is the purpose of this paper to solve generalised strictly hyperbolic problems, which might 
be generated by singular coefficients and data, by means of FlO-techniques placed in the Colombeau 
context. Generalised pseudodifferential and Fourier integral operators as well as the microlocal tools in 
the Colombeau framework that have been developed over the past years, provide new and powerful tools 
and techniques which were not available at the time of the first work on hyperbolic generalised systems 
and equations in [18, 22, 25, 26]. On the other hand, the application of FIOs to generalised strictly 
hyperbolic problems is very natural also in the Colombeau context and will have interesting connections 
with systems with multiple characteristics and singularities. This and the singularity structure in the 
case of systems and of higher order equations is the subject of ongoing research and will be published 
elsewhere. 

We comment on some technical aspects of the theory. Elements of the Colombeau algebra C/(IR'"+^) are 
equivalence classes of nets {ue)e of C°^-functions satisfying asymptotic bounds of order in tcirms of 
the local L°°-seminorms of all derivatives as £ — > 0. The generalised wave front set is also defined by 
means of decay in the dual variable and uniform asymptotic bounds as £ — > 0. The existence theory is 
based on energy estimates and Gronwall's inequality and requires more restrictive asymptotic bounds of 
the type |log£] at certain places; propagation of regularity may require additional so-called slow scale 
estimates. The necessity of such stronger bounds has been argued in [11, 19, 26]. When the coefficients 
are discontinuous functions or distributions, corresponding members of the Colombeau algebra satisfying 
the stronger bounds can always be constructed by suitable regularisation [25]. Admitting distributional 
coefficients in hyperbolic equations may lead to infinite propagation speed [22] . To avoid this phenomenon, 
one may assume that the coefficients are constant for large \x\, as is often done in the classical case as 
well. Matters are even more simplified by assuming that the coefficients are compactly supported in x. 
This still allows one to model and study singularities as strong as desired (in the non-constant regime) , 
but facilitates the application of Sobolev estimates. 

It may also be useful to employ Colombeau spaces with asymptotic estimates based on Sobolev spaces. 
It has been shown in [18] that the first order scalar hyperbolic pseudodifferential problem 

(0.2) dtU-ia{t,x,Dx)u = f, u{0) = g, 
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with a real valued principal part has a unique solution in the Colombeau algebra t/2,2([0,T] x R") based 
on H°° . Assuming that the symbol a is compactly supported in a;, the Sobolev embedding theorem yields 
existence and uniqueness of a Colombeau solution in S([0,T] x R") as well. 

Proofs of the mentioned existence and imiquencss results are based on a priori L^-estimates in the 
representing nets of smooth functions. These methods have been successful regarding the well-posedness 
of the Cauchy problems under consideration but do not provide deep qualitative information on the 
solution. This already applies to the classical case of hyperbolic systems with C^-coeiJicients and has 
been one of the main motivations for the constructive FIO approach introduced by Duistermaat and 
Hormander [5, 6, 17], which we extend here to the Colombeau framework. For a survey on local and 
global regularity of Fourier integral operators on L'p spaces as well as on Colombeau spaces we refer the 
reader to [29]. 

The plan of the paper is as follows. In Section 1, we recall required notions from the Colombeau theory 
of generalised functions, in particular, microlocal tools, generalised symbols and phase functions and 
generalised FIOs. Section 2 is devoted to scalar first order hyperbolic partial differential equations 
with Colombeau coefficients. We show how the generalised solution can be represented by means of a 
generalised FIO. This includes the construction of Colombeau solutions to the eikonal equations and to 
the transport equation for the generalised symbol. In Section 3, we treat the case of first order hyperbolic 
pseudodifferential equations whose principal part depends on time only. This has the advantage that 
generalised solutions to the eikonal equation can be given explicitly (the more intricate t- and cc-dependent 
case is postponed to work in preparation on hyperbolic problems with singularities). We solve an infinite 
system of transport equations yielding the asymptotic expansion of the generalised symbol, and construct 
an FIO parametrix of the gcineralised solution. Section 4 addresses microlocal properties of the kernel of 
a generalised FIO as well as its action on a Colombeau generalised functions. This allows one to compute 
the generalised wave front set (in space and time) of solutions to first order hyperbolic partial differential 
equations with non-smooth coefficients depcniding on time. Further, the spatial wave front set of solutions 
to first order partial differential equations at fixed time is computed, when the coefScients are generalised 
functions depending on space and time, in case the Hamiltonian ffow has a limit. We show that the 
generalised wave front set is invariant under the Hamiltonian flow for time-dependent coefRcients. The 
paper concludes with some explicit examples involving jump discontinuities and delta functions in the 
coefficients. 



1 Basic notions 

1.1 Basic notions of Colombeau theory 

This section gives some background on Colombeau techniques used in the sequel of this paper. As main 
sources we refer to [8, 11, 12, 15]. 

Nets of complex numbers. A net {ue)e in C^"'^! is said to be strictly nonzero if there exist r > 
and 77 € (0, 1] such that ]ug] > e"^ for all e G (0,77]. For several regularity issues we will make use of the 
concept of slow scale net (s.s.n). A slow scale net is a net {rg)g e C'"'^! such that 

Vg>0 ]r^]9 = 0(e-i) as e 0. 

A net (?ie)e in C'^"'^' is said to be slow scale- strictly nonzero if there exist a slow scale net {ss)s and 
7? e (0, 1] such that \ue\ > l/sg for aU e G (0,??]. 

C-modules of generalised functions based on a locally convex topological vector space. The 

most common algebras of generalised functions of Colombeau type as well as the spaces of generalised 
symbols we deal with are introduced by referring to the following general models. 

Let E he a locally convex topological vector space topologised through the family of seminorms {pi}iei- 
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Vie/3iVeN pj(ue) = 0(e"^)ase ^ 0}, 
Vi e / 3{uje)e s.s.n. Pi(ue) = ©(cje) as £ -> 0}, 
3N€N\/iGl p,{u,) ^ O {e-'^) as e^O}, 
VieJVgeN = 0(e«)ase^O}, 

are called _E-nioderate, E'-moderate of slow seale type, i?-regular and i^-negligible, respectively. We define 

the space of generalised functions based on E as the factor space Qe Me/J^e- 

The ring of complex generalised nvmhers, denoted by C, is obtained by taking E = C Note that C is not a 
field since by Theorem 1.2.38 in [15] only the elements which are strictly nonzero (i.e. the elements which 
have a representative strictly nonzero) are invertible and vice versa. Note that all the representatives of 
u G C are strictly nonzero once we know that there exists at least one which is strictly nonzero. When u 
has a representative which is slow scale-strictly nonzero we say that it is slow scale-invertible. 

For any locally convex topological vector space E the space Qe has the structure of a C-module. The C- 
module := M^eI^e of slow scale regular generalised functions and the C-modulc := M'^ /Me of 
regular generalised functions are subspaces of whose role is to describe different notions of regularity. 
We use the notation u = [{ue)e\ for the class u of (we)e in Qe- This is the usual way adopted in the paper 
to denote an equivalence class. 

The most common Colombeau algebras. The Colombeau algebra Q{0) = £m{^)/M{^) (see [15]) 
can be obtained as a C-modulc of ^_E-type by choosing E = C°°{il). From a structural point of view 
il — > Q{fl) is a fine sheaf of differential algebras on R". ^c(^) denotes the Colombeau algebra of 
generalised functions with compact support. 

Regularity theory in the Colombeau context as initiated in [26] is based on the subalgebra Q°°{Sl) of all 
elements u of Q{il,) having a representative (me)^ belonging to the set 

£^{n) := {{ue)e e £[n] : VK (sClBNe NVa e N" sup \d"u^{x)\ = 0(£-^)as e 0}. 

xeK 

The 0°°-singular support of w G Q(fl) (singsuppgoc u) is defined as the complement of the set of points 
X such that u\v € Q°°{V) for some open neighborhood V of x. 

The intersection Q°°{Sl) Ci Qd^) will be denoted by Q^{Q,). In the course of the paper, for issues related 
to the Fourier transform, we will make use of the Colombeau algebra ^(R") = Qy(^n) of generalised 
functions based on ^(R") and of the corresponding regular version Q^{M."). Finally, we recall that the 
Colombeau algebra ^^(IR") of tempered generalised functions is defined as fT-(R")/A/'r(K"), where fT-(IR") 
is the space 

{{us)e G Om(K")(°'11 : Va G N" 3iV e N sup (1 + |a;|)-^|a«Ue(a;)| = Cl(e"^) as e 0} 

of r-moderate nets and ^/^(R") is the space 

{{ue)e G Om(K")(°'^1 : Va G N" 3Ar e NVg G N sup (1 + \x\)-^\d'^udx)\ = 0(£«) as e 0} 

of r-negligible nets. The subalgebra Q^{M.") of regular and tempered generalised functions is the quotient 
£:^(K")/A/;(R"), where £"~(R") is the set of all {ue)e G 0m(K")'°'^1 satisfying the following condition: 

3N G NVa G N" 3M G N sup (1 + |x|)-^|5«Ue(a;)| = 0(e-^). 

xeR" 

Special types of Colombeau generalised functions in ^(R"). An element of the Colombeau algebra 
is real valued if it is represented by a net of real valued smooth functions. We recall that u G ^(R") is 



The elements of 

Me := {(Me)e G : 

:= {(u,), e : 

Me := {(ue)e G : 
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- of logarithmic type if there exists a representative {ug)g with the property that for ah K (g R" 

sup |w£(a;)| = 0(log(£~^)) as £ — >■ 0; 

xeK 

- of slow scale logarithmic type if there exists a representative {ue)s with the property that for all 

g M" there exist a slow scale net {He)s such that 

sup \us{x)\ = 0(log(/ie)) as e ^ 0; 

- slow scale regular if there exists a representative {us)s with the property that for all K (g R", for 
all a e N" there exist a slow scale net (/ie)^ such that 

sup \d°'ue{x)\ = 0{iie) as £ ->■ 0. 

xeK 

Clearly the previous properties hold for all representatives of u once they are known to be valid for one. 

Topological dual of a Colombeau algebra. A topological theory of Colombeau algebras has been 
developed in [7, 8]. The duality theory for C-modules presented in [8] in the framework of topological and 
locally convex topological C-modules, has provided the theoretical tools for dealing with the topological 
duals of the Colombeau algebras Gc{^) and We recall that £{Q{fl),C) and £(^c(^^))C) denote 

the space of all C-linear and continuous functionals on Q{il) and Qc{^), respectively. For the choice of 
topologies given in [7] one has the following chains of continuous embeddings; 

(1.3) g°°{Q) c g{n) c £{g^{n), c), 

(1.4) g^m c gd^) c c{g{n),c), 

(1.5) £(e?(n),c) c£(e?,(f2),c). 

In (1.3) and (1.4) the inclusion in the dual is given via integration (u — !> (v — !> J^-^u{x)v{x)dxj'j (for defi- 
nitions and properties of the integral of a Colombeau generalised functions see [15]) while the embedding 
in (1.5) is determined by the inclusion Gd^) CI g{^). Since fl — > £(^7c(f^),C) is a sheaf we can define 
the support of a functional T (denoted by supp T). In analogy with distribution theory we have that 
C{g{n),C) can be identified with the set of functionals in C{gc{^),C) having compact support. 

For qiiestions related to regularity theory and microlocal analysis particular attention is given to those 
functionals in £(C/c(f^)jC) and C{g{il),C) which have a "basic" structure. In detail, we say that T e 
C{gd'^),C) is basic if there exists a net {T^)e € D'(0)(°'il fulfilling the following condition: for all K Q. 
there exist j €N, c> 0, N gN and ?? G (0, 1] such that 

V/Gl?K(O)V£e(0,?7] \Tdf)\<ce-^ sup \d"f{x)\ 

xeK,\a\<j 

and Tu = [{T^Ue)e] ioi alluegdi^). 

In the same way a functional T G jC{g{fl),C) is said to be basic if there exists a net (Te)^ G £'(0)'°'^! 
such that there exist K(S:fl, j£f^,c>0, Ne.N and rj G (0, 1] with the property 

V/GC°°(n)V£G(0,77] |Te(/)|<c£-^ sup \d"f{x)\ 

and Tu = [{TsUs)e] for all u G g{fi). The sets of basic functionals on gd^) and S(ri) are denoted by 
£b(^c(fi),C) and Ch{g{'^),Q, respectively. 

The Fourier transform on ^(M"), /:(^(M"), C) and jC{g{n),C). The Fourier transform on ^(R") 
is defined by the corresponding transformation at the level of representatives, as follows: 

J- : ^(R") ^ ^(R") : u ^ [{uiU 
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is a C- linear continuous map from ^(R") into itself which extends to the dual in a natural way. In 
detail, we define the Fourier transform of T e £(^(]R"),C) as the functional in £(^(M"),C) given by 

J^{T){u)=T{Tu). 

As shown in [9, Remark 1.5] C{Q{ft),C) is embedded in £(g^(M"),C) by means of the map 
£{g{n), C) ^ £(e^(M"), C):T^{u^ T{{u,\n)e + m^))) ■ 

In particular, when T is a basic functional in C(g(fl),C) wc have from [9, Proposition 1.6, Remark 1.7] 
that the Fourier transform of T is the tempered generalised function obtained as the action of T[y) on 
e-*f«, i.e., J"(r) = T{c-''<) = {T^{c-'<))e + A/;(R"). More precisely, J"(T) belongs to g^{R"). 

Microlocal analysis in the Colombeau context: the ^°°-wave front set for generalised func- 
tions and functionals. For an introduction to microlocal analysis in the Colombeau context we refer 
to [9, 12]. Here we only recall those microlocal concepts which we will employ in the final section of the 
paper. The ^°°-wave front set of u € Qd^) (WF'gu) is defined as the complement of the set of points 
(.To,^o) G T*{ft) \0 fulfilling the following property: there exists a representative (u^)^ of u, a cut-off 
function e C^{^) with (f{xo) = 1, a conic neighborhood T of and a number N such that for all 
I e R 

sup = 0{e-^) as £ 0. 

By construction TrnWFgooU = singsuppgoo u. In addition, Theorem 3.11 in [12] shows that WFgoo 
coincides with the set Wci,e°o(u) := C\Aueg°°(Q)^^^^(^)' where the intersection is taken over all the 
standard properly supported pseudodifferential operators A G \I'°(f2) such that Au € Q°°{^1). The 
adjective standard refers to symbols which do not depend on the parameter s but belong to the usual 
Hormander classes. This pseiidodifferential characterisation provides the blueprint for extending the 
notion of wave front set from g{Q) to the dual £(^c(^),C). In detail, for T G £(5c(^),C) we define 

(1.6) WFgoo(T) = Wci,5oc(T) := f] Char(yl) 

ATeg°°{n) 

where the intersection is taken over all the standard properly supported pseudodifferential operators 
A e such that AT e g°°{Q.). It follows that 7ra(WFgooT) = singsuppgoo T. 

A useful characterisation of WFgoo (T), in terms of estimates at the Fourier transform level, is valid when 

T is basic. It involves the sets of generalised functions ^^Q(r), where F is a conic subset of R" \ 0, of all 
tempered generalised functions u having a representative (u^)^ fulfilling the following condition: 

3NeN\/leR sup{0'\ue{0\=O{e-^) as £ -> 0. 

Let T e C{Qc{^), C). Theorem 3.13 in [9] shows that {xo,^o) ^ WFgooT if and only if there exists a conic 
neighborhood F of and a cut-off function ip e with ip{xo) = 1 such that T{'^T) e G'^q(X)- 

1.2 Generalised Fourier integral operators: generalised symbols and phase 
functions 

In this subsection we collect some basic notions concerning generalised pseudodifferential and Fourier 
integral operators. For a detailed presentation we refer to [10, 13]. 

Generalised symbols. Let be an open subset of M". By a generalised symbols of order m we mean 
an element of the space Ge based on E = S'^^i^ x W). Analogously, with E = S'^giCt x Rp) is the 
space of slow scale regular generalised symbol. We say that a generalised symbol a of order m is regular 
if it has a representative {a^)^ fulfilling the following condition: 

yK ^flBN emj eN |oe|^] = sup |a,|^^ . = 0(£-^) as£^0. 

H+W\<j 
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A notion of asymptotic expansion for generalised symbols has been introduced in [10, Subsection 2.5] 
and employed in developing a complete symbolic calculus for generalised pseudodifferential operators, see 
[10, 11]. Finally, we recall that the conic support of a generalised symbol a of order m is the complement 
of the set of points (xq, ^o) G ri x R^* such that there exists a relatively compact open neighborhood U of 
xo, a conic open neighborhood F of and a representative (ae)^ of a satisfying the condition 

(1.7) VaeNfV/3GN"V9eN sup {0~'^'^^^°'^~^^^Wd^a,{x,0\ = Oie") as e ^ 0. 

xeu,ier 

By definition, conesuppa is a closed conic subset of O x R^. 

Generalised phase functions. A phase function (l){y, ^) on x R^ is a smooth function on x R^ \ 0, 
real valued, positively homogeneous of degree 1 in ^ such that Vy^((p{y, 7^ for ally Gfl and ^ e K^\0. 
We denote the set of all phase functions on x R^ by $(0 x W) and the set of all nets in $(0 x Rp)*"'^! 
by $[r2 X W]. We recall that Sl^{il x \ 0) is the space of symbols on x R^ \ homogeneous of order 
1 in ^, i.e. sup^g;f^jgRp\o \^\~^^"\9^dxaix,0\ < oo for all K Q, for all a G and /3 e N". 



Definition 1.1. An element of A4<s>{Cl x R^") is a net 

(i) {(j)e)e e Msl^(ny.«.v\0), 

(a) for all K the net 



e $[il X 



satisfying the conditions: 



inf 



1^1 



is strictly nonzero. 



On Ai^{ClxMP) we introduce the equivalence relation ~ as follows: {(t>e)e 
Afgi (nxMP\o)- elements of the factor space 

$(n X W) := M^{n X 
will be called generalised phase functions. 



{uje)e if and only if {(j)e-uje) G 



Finally, we say that 4> G x R^*) is a slow scale generalised phase function if it has a representative 
((/>e)e G Ad^gi (QxMP\o) ^'^^^ ^^^^ (*^) scale strictly nonzero. 



In the sequel W is an open subset of R" . We denote by Q x R^] the set of all nets {4'e)ee{o,i] of 
continuous functions on fi' x O x R*" which are smooth on O' x O x R^' \ {0} and such that {(j^eix, •, •))£ € 
$[0 x W] for all x G f2'. 

Definition 1.2. An element o/A^$(0';0 x R^') is a net {(j>s)e G x R^'] satisfying the conditions: 

(i) {4'e)e € Msl^(n'y.nxm.f\0)> 

(a) for all K' (E O' and K <s^Cl the net 



(1.8) 

is strictly nonzero. 



inf 



On X W) we introduce the equivalence relation ~ as follows: {<pe)e ~ ('^e)e */ '^"'^ only if 

((/>e — iJs)s € A/gi (o'xOxRpXo)- ^'^fi elements of the factor space 

^n';n X RP) := M^{n';n X RP)/ - . 

are called generalised phase functions with respect to the variables in f2 x R^. If (p € $(il'; x R*") has 
a representative {(j)e)e G (n'xOxRp\o) ^^^'^ ^^^^ '"'^^ ^'^ i^-^) ■^'^^ scale strictly nonzero then it 
is called slow scale generalised phase function with respect to the variables in x R^. 
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Slow scale critical points. 

Definition 1.3. Let (j) G $(0 x Rp). We define C ^ x W \ as the complement of the set of all 
{xo,^o) € fl X MP \ with the property that there exist a relatively compact open neighborhood U{xo) of 
Xq and a conic open neighborhood r(^o) C MP \ o/ such that the generalised function |V|(?i(-, •)|^ is 
slow scale-invertible on U{xq) x r(^o)- We set iraiC^) = and R^^ = {S^^y. 

By construction is a conic closed subset of x IRp \ and R^^ C R^ C is open. It is routine to 
check that the region coincides with the classical one {(x,^) S x M'' \ : V^(l){x,£,) = 0} when </> 
is a standard phase function independent of e. 

Generalised Fourier integral operators. Let 4> G x MP), a G Gs'"^{s^'xQxrp) and u G Qc{^)- 

The generalised oscillatory integral 

I^{a){u){x) = [ e'^(^'^'«a(:r, y, Ou{y) dy 

defines a generalised function in Q{fl') and the map 

(1.9) A : ^?e(0) ^ g(O') : u ^ /^(a)(«) 

is continuous. The operator A defined in (1.9) is called generalised Fourier integral operator with am- 
plitude a and phase function (p G ^{ft';fl x M^). From Theorem 4.6 in [10] wc have that when phase 
function and amplitude are both slow scale regular the corresponding generalised Fourier integral oper- 
ator maps ^^(O) continuously into Q°°{Q'). If a e Gs'-°°{n'xnxRp) ^^^^ ^ maps Q^i^) into ^°°(0'). 
Pseudodifferential operators are special Fourier integral operators with ^{x,y,^) = {x — y)^ in (1.9). 

Composition of a generalised Fourier integral operator with a generalised pseudodifferential 

operator. We conclude this survey on generalised Fourier integral operators by studying the composition 
with a generalised pseudodifferential operator. First of all we focus on Fourier integral operators of the 
form 



FUb){u){x)= [ e'"^-'^^b{x,v)uiv)'^V, 



where w e QsI^{Q'xR''\o)^ b G ^s"'(a'xM") and u G Qc{^). Note that (j){x,y,r]) := cv{x,r)) - j/r? is a well- 
defined generalised phase function belonging to $(0';n x M") and 7,^(6) = -FL(6). Theorem 5.11 in [10] 

provides the composition formula stated below. 

Composition formula for a{x, D)F^{b). Let oj G G^gi (qxM"\o) have a representative {oje)^ such that 
V^oJe ^ for all e G (0, 1] and for all K Q 



inf 



is slow scale strictly non-zero. Let a G 55m(QxK'') ^ ^ ^s'(nxK''\o) ^^^^ supp^ b <s il. Then the 
operator a{x,D)F^j{b) has the following properties: 

(i) it maps g^{n) into g°°{n); 

(ii) it is of the form 

e"^(^'^)/i(a;, ri)u{ri) drj + r{x, D)u, 



I 



where h G ^5i+m(nxM"\o) asymptotic expansion given by the symbols 

d?a(x,VxOJix,r])) , , 
K{x, v) = ^ J D'^ (e-(^'-''')6(2, ry)) a G N", 

with uj{z,x,r]) := uj{z, ri)—uj{x, rf) — (Vxw(a;, ri),z — x), and r is slow scale regular and of order — oo. 
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2 Transport equations with generalised coefficients 



In this section, we are concerned with the Cauchy problem for the first order hyperbohc equation 

n 

-DtM = > tti i(t,x)DiU + ar)(t,x)u 

(2.10) fr[ 

u(0, •) =uo, 

where Dj = D^^, the coefficients aij are real valued Colombeau generalised functions in ^(R"+^), 
do & ^(M""'"^) and uq S tJc(IR")- The following theorem combines the well-posedness results of [22] with 
the more recent investigations of [18]. 

Theorem 2.1. Let the coefficients aij, j = 1, n and uq be Colombeau generalised functions in ^^(R""''^) 
compactly supported in x. Assume that the coefficients aij are real valued and dkaij as well as ao are 
of logarithmic type (k,j=l,...,n). Then: 

(i) For each uq S ^(M") the Cauchy problem (2.10) has a unique solution u S ^(R"+^). 
(ii) If in addition uq € Sc(K") then the solution u is compactly supported in x. 

(Hi) If in addition ao as well as uq are real valued then the solution u is a real valued generalised function. 

(iv) If the coefficients ai_j and oq are slow scale regular and dkai j (k,j=l,...,n) as well as oq are of slow 
scale logarithmic type then for each uq G ^°°(M") the unique solution u € ^(M"+^) of the Cauchy 
problem (2.10) belongs to g°°{W+^). 

(v) Under the hypotheses of {iv) if the initial data uq are slow scale regular then u is slow scale regular 
as well. 

The aim of this section is to prove that the solution u of the Cauchy problem (2.10) can be written as 
the action of a generalised Fourier integral operator F^{b) on the initial data uq € Qd^'^)- This requires 
to determine the phase function (p and the symbol b. 

2.1 The generalised phase function and the characteristic curves 

The generalised phase function ^ is the solution of the eikonal equation determined by the principal part 
of the operator 

n 

Dt - '^aij{t,x)Dj - ao{t,x) 
under the initial condition (j){Q, x, rf) = xr}. Thus one has to solve the linear Cauchy problem 

n 

dt<i){t,x,ri) =^ai^j{t,x)dj<j){t,x,ri), 
(2-11) j=i 

(/)(0, X, ?]) = xrj. 

Under the hypotheses of Theorem 2.1 on the coefficients a\j, we already know that there exists a unique 
solution (f) G Q{W^^'^). More precisely, it has the form 

n 

<i){t,x,ri) = ^ujh{t,x)rih, 
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where u;^, = 1, n, are solutions of the Cauchy problems 

n 

(2-12) U 

In the following proposition we describe the properties of the generalised functions w/j more specifically. 
The solutions 71, ...,7„ of the initial value problem 

(2.13) ^7h(x,i,s) = -ai,,j(s,7i(a;,t,s),72(x,t,s),...,7„(x,t,s)), 'yh{x,t,t) = Xj, h==l,...,n. 

are the components of the characteristic curve 7 = (71, ...,7„) associated with the differential operator 
^J^j^ fli J (t, x)_Dj. Note that from Theorem 1.5.2 and Remark 1.5.3 in [15] this initial value problem 
is well-posed in when aij is compactly supported in x and has first order cc-derivatives of 

logarithmic type. 

We make the following assumptions: 

{hi) the coefficients aij are real valued generalised functions in 5(]R"+^), compactly supported with 
respect to x with dkaij of logarithmic type {k,j = 1, ...,n); 

(/i2) the coefficients aij are real valued slow scale regular generalised functions in compactly 
supported with respect to x with dkaij of slow scale logarithmic type (k,j=l,...,n). 

The compact support property is introduced here to keep the presentation simple and could be relaxed. 
Proposition 2.2. 

(i) Under the hypothesis (hi) there exists a unique real valued solution ujhit^x) £ QiMP'^^) of the 
Cauchy problem (2.12); ujh{t,x) is the h-th component of the characteristic curve ^{x,t,0). 

(a) Under the hypothesis {h2) the solution ujh is slow scale regular. 

Proof. From Theorem 2.1 (i) and (Hi) it is clear that there exists a unique real valued Colombeau 
solution ojh € ^(M"+^). Since the initial data x^ are smooth and therefore slow scale regular we have 
from Theorem 2.1 {iv) and (v) that uih is slow scale regular under the set of assumptions {h2) on aij 
and dkUij, k,j = 1, ...,n. It remains to prove that uJh{t,x) = 7?i(a;,t,0). This comes from the fact that 
coh is constant along the characteristic curves ^{x,t,s), i.e., working at the level of representatives, 

-^i^hA^^ 7i,e(a;, t, s), 72,e(a;, t, s), 7„,e(a;, t, $)) = 0. 

Hence, 

i^h,e{t, 7i,s{x, t, t), 7„,£(a;, t, t)) = ujh,e{0, 7i,e{^, t, 0), 7„,£(a;, t, 0)) 
for each t G R. This implies that ujhit, x) = jh{x, t, 0) in ^(M"+i). □ 

Summarising we can state the following proposition. 

Proposition 2.3. // {hi) holds then the generalised phase function 

n n 

<t>{t,x,v) =^ujh{t,x)rih = ^'yh{x,t,0)r]h 

h=l h=l 

solves the eikonal Cauchy problem (2.11). 

If {h2) holds then (j) is a slow scale generalised phase function. 
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Remark 2.4. The Colombcau generalised function is actually a generalised symbol homogeneous of 
order 1 in t/. With a certain abuse of language we employ the expression generalised phase function, 
previously referred to <p' {t,x,y,rf) = 4>{t,x,rf) — yrj. Indeed, for 0' one has the typical invertibility 
condition on the gradient, i.e., 

\V(t)'{t,x,y,r])\ = |(Vt^(f,a;,77), V,:<^(i,a;,77),-?7, V,,(/>(t,a;,r?) -y)\>l 

for all 7] with |r7| = 1. 

2.2 The transport equation for the generahsed symbol 

To compute the generalised symbol b we need to solve the Cauchy problem (2.10) with initial condition 
1 at f = 0. In detail, 

n 

Dtb = oi -j(t, x)Dib + aoit, x)b 
(2.14) ^ 

6(0, •) = 1. 

We introduce the following set of hypotheses on ao- 

(11) ao is a generalised function in ^(R"+^), compactly supported in x and of logarithmic type; 

(12) ao is a slow scale regular generalised function in ^(1R"+^), compactly supported in x with 0-derivative 
of slow scale logarithmic type. 

Using Theorem 2.1 and integrating along the characteristics we have the following existence and unique- 
ness result. 

Theorem 2.5. 

(i) Under the hypothesis (hi) on the coefficients aij and the hypothesis {il) on oq there exists a unique 
solution b € 5(M"+i) of the Cauchy problem (2.14). 

(a) Under the hypothesis (/i2) on the coefficients aij and the hypothesis {i2) on ao, the solution b G 
^(IR"'+^) of the Cauchy problem (2.14) is slow scale regular. 

(Hi) Further, b{t,x) = e^'^^*'^) with 

/3{t,x)= / ao{s,'yi{x,t,s),...,j„{x,t,s))ds. 
Jo 

2.3 Generalised FIO formula 

A combination of Proposition 2.3 with Theorem 2.5 yields the following FIO formula. 
Proposition 2.6. 

(i) Under the hypotheses (hi) and (il) on the coefficients aij and ao, respectively, the solution u G 
QiW^^^) of the Cauchy problem, (2.10) can be written as 

(2.15) u{t, x) = F^{b)(uo){t, x) := [ e^*(*'^''')6(t, x)u5{rj) drj, 

where the generalised phase function cj) is defined by Proposition 2.3 and b £ ^(]R"+^) by Theorem 
2.5. 

(ii) Under the hypotheses (/i2) and {i2) the formula (2.15) holds with (j) and b slow scale regular. 
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2.4 The non- homogeneous Cauchy problem 

We conclude this section by finding a solution formula for the non-homogeneous Cauchy problem 

n 

DfU = y^ai j{t,x)DjU + ao{t,x)u + f{t,x) 
(2.16) ' 

u(0, •) =uo, 

where / G tJ(K""'"^) is compactly supported with respect to x. Note that the Fourier integral operator 
F^{b) of Proposition 2.6 is given by 

F^{b){uo){t,x) = b{t,x)uo{j{x,t,0)). 

It defines, for each t e M, a map 

Uit) = F^{b){t) : GciW) ^ ^c(R") : uo ^ F^{b){uo){t, •) 

such that U{Q) = I and 
Theorem 2.7. 

(i) Under the hypotheses (hi) and (il) on the coefficients aij and ao, respectively, the solution u G 
of the Cauchy problem (2.16) can be written as 

(2.17) u{t,x)=F^{b){t)(uo + iJ^ ^(^^^(^^o^^)) /(^.7(-,0,r))c;r^(a:) 

where the generalised phase function (j) is defined by Proposition 2.3 and b G ^(IR"+^) by Theorem 
2.5. 

(a) Under the hypotheses {h2) and (i2) the formula (2.17) holds with <j) and b slow scale regular. 

3 First order hyperbolic pseudodifferential equations with ge- 
neralized symbols 

We now consider hyperbolic first order pseudodifferential equations of the type 
(3.18) Dtu ^ ai{t,Dx)u + ao{t,x,Dx)u 

where ai and oq are generalised symbols of order 1 and 0, respectively, with ai real valued and independent 
of X. As mentioned in the Introduction, we restrict ourselves to ^-dependent principal parts. We begin 
by collecting what is known about equation (3.18) in the Colombeau context. The following theorem is 
due to Hermann in [18]. The well-posedeness of the Cauchy problem 

Dtu = ai{t,D^)u + ao{t,x,D^)u + f 
u{0,-)=Uo, 

is intended in the Colombeau algebra 5_f/°°((-T,T)xE") based on H°°{{—T,T) x K"). Here we use the 
notation ^2,2((^r, T) xM") introduced in [2]. The choice of this setting is motivated by a uniqueness issue: 
the solution u to the problem (3.19) fails to be unique in the usual Colombeau algebra g{[—T,T] x R") 
whereas is uniquely determined in ^2,2 ((—?", T) x R"). Finally, with the expressions generalised symbol 
and slow scale regular generalised symbol we refer to the elements of the spaces ^c=°([-t.t],S"»(r2")) and 
^c'="([-TT] S"*(R2"))' respectively. Note that with the notation S"^{M?") we intend symbols satisfying 
uniform estimates with respect to x G M" and ^ G M", i.e., sup(^_j)gR2,i(0~'"+'"'|9|'c^f a(a;,^)| < 00. In 
the sequel {k, h) denotes the ?7-derivatives and a;-derivatives of a symbol up to order k and h, respectively. 
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Theorem 3.1. Let ai be a real valued generalised symbol of order 1, uq a generalised symbol of order 0, 
/ e a2,2((-T,T) X M") and uo e a2,2(M"). 

(i) If there exist representatives {ai^e)e md (ao,e)£ of ai and gq respectively such that {ai^e)e is of 
log-type up to order + with A:„ = 3([n/2j + 1) anc? Z„ = 2(n + 2), and (ao.e)^ is of log-type 
up to order {k[^,jn) with k'^ = l'^ = [n/2j + 1, then the Cauchy problem (3.19) is well-posed in 

02,2 ((-T, T) xM"). 

(ii) If for large \x\ the net {ao^e)e does not depend on t then one can set fc„ = 1, l„ = n + 2, fc^ = 0, 
l'^ = n-\-l in {i). 

(Hi) If ai and ao are slow scale regular generalised symbols and the log-type conditions in (i) on {ai^e)e 
and (ao,£)e are replaced by slow scale log-type assum,ptions then for each f G CJ|°2((~2^) ^) x K") and 
uq e ^2^2(^^"') unique solution u to the Cauchy problem (3.19) belongs to Q2^2{{~T,T) x M"). 

Our aim is now to find an FIO formula by which to express the solution u. More precisely, we will 
construct a generalised FIO parametrix for the hyperbolic Cauchy problem 

Dtu = ai{t,Dx)u + ao{t,x,Dx)u 

^ ■ ' u(0,-) = woegc(M"). 

Under suitable moderateness assumptions we will get well-posedness and -regularity. 



3.1 The generalised phase function and the eikonal equation 

Wc begin by determining the generalised phase function (l){t,x,r]), i.e. by solving the following eikonal 

equation. 

Proposition 3.2. The eikonal equation 

^2 dt(t){t,x,r]) = ai{t,V^(t){t,x,r])), 

(p{0,x,r]) = XT], 

has the solution 

I ds 



(p{t,x,r]) =xr]+ ai{s,r]) 
Jo 



(i) TO ^C~([-T,T], Si(K^")) ^/ «1 € ^C~([-T,T], Si(R2")), 

(ii) in G^C'=°([-T,T\,s^(VL^^)) ^■^'^i ^ ^c~([-t,t],si(r2"))- 

Remark 3.3. (i) Note that (j){t, x, rj) is a symbol of order 1 with Va;</'(i, x,ri) = rj but non-homogeneous 
with respect to r]. This is not an obstacle in defining the Fourier integral operator 

f ei'l>(*'^''n)h{t,x,r])v{r])dri= [ e''-'^^*''''''^-y''^b{t, x,T])viy) dy d^rj 

for V G 0c (R") and a generalised symbol b. Indeed, for giving a meaning to the oscillatory integral on 
the right-hand side (by means of an operator to be used in the integration by parts), it is sufficient 
that the following bound from below 

(3.22) \V,<l>',{t,x,r])f + \Vycl>',{t,x,r,)f + \rjf\\7^cf,',{t,x,v)f > XM^ 

holds for all t G [— T, T], {x, y, t]) G M?^ and e G (0, 1], with 0'(<, x, y, rj) = 0(i, x, rj) — yq and some strictly 
non-zero net Ag. The condition (3.22) is trivially satisfied by the phase function in Proposition (3.2). 
Indeed, for 



(f)' {t, x,y,ri) = {x-y)r]+ I ai (s, rj) ds 

Jo 
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one has |Va;0^(t, x, = jryp. 

(ii) Easy computations at the level of representatives show that the Fourier integral operator with phase 
function (/>(t, x, 77) = xr] + a-i{s,r]) ds and symbol b{t,x,r]) £ Qc'^{[-T,T],S"^{m?'^)) maps ^c(K") (or 
^(R")) into e'c<»([-T,T],j^(M"))- Recalling the embedding g^iRn) C ^(R") C ^2,2 (R") (see [8]) we have 

F^(6) : ^(R") ^ e?coc([-T,T],^(K")) ^ e?2,2((-T,T) x R") 

and 

ai{t,D^)F^{b) + ao{t,x,D^)F^{b) : a^(R") ^ acoc([_T,T],^(Mn)) ^ a2,2((-T,T) x M"). 

In addition, if we work with slow scale generalised phase functions and symbols then the previous mapping 
properties hold between e~(R") and g^2{{-T,T) x R"). 

Theorem 5.10 in [10] concerning the composition of a generalised pseudodifferential operator with a 
generalised Fourier integral operator can be easily adapted to Fourier integral operators with a phase 
function as above. This means that under the slow scale assumptions of Theorem for any 

b G ^c''°°([-T,T],S"'(R2n)) and V e ^c(]R"), we have 

(3.23) ai{t,D^)F^{b)v + ao{t,x,D^)F4b)v= [ e''''^*-'=-'''^h{t,x,r])v{r]) drj + r{t,x, D^)v, 

where h and r arc generalised symbols of order m + l and —00, respectively, v G ^(R") and the equality 
is intended in ^2,2 ((—r,T) x R"). If we take slow scale regular phase functions and symbols then h and 
r are slow scale regular as well and r{t,x,Dx)v G g22{{—T,T) x R"). The symbol h has the following 
asymptotic expansion (as defined in [10]): 

MM,.)~ E ^"(*.-.'^)= E ^^"(-^(^'-^) + -»(^'-'-^)) ^g(e-^6(M..))U.. 



= E ^^°("^(^'^)+^°(^'^''^)) z)g6(M,ry). 
3.2 The transport equations and the generaUsed symbol 

We proceed in the construction of a generalised FIO parametrix for the Cauchy problem (3.20) by looking 

forasymbolfo G ^c~([-T,r],so(M2")) given by the asymptotic expansion X^feeN ^fe' G ^^c~([-T,T],s-'i(M2'»))- 
The symbols bk will solve some specific transport equations and will determine b such that 

DtF^{b)-ai{t,D^)F^{b)-ao{t,x,D^)F4,{b) 

is a ^2°2"'"6gula'rising operator. We begin by observing that by the composition formula (3.23), 

ai{t, D,)F^{h)v + ao(t, x, D,)F^{b)v = F^{h)u + r{t, x, D,)v, v G Gd^"), 

where r has order —00 and is slow scale regular when ai and ao are slow scale regular as well. In other 
words we can write 

DtF4b)v - ai{t, D^)F^{b)v - ao{t,x, D^)F^{b)v = DtF^{b)v - F^{h)v - r{t, x, D^)v. 

Note that taking a suitable cut-off function tp we can write the action of r{t,x,Dx) on v as the integral 
operator 

xl}{y)kr{t,x,y)v{y)dy, 



j 



with ip{y)kr{t, X, y) G S2°2((~^; T) xR^") rapidly decreasing with respect to x and y. Then the generalised 
function j^2n'^{y)kr{t,x,y)v{y) dy belongs to ^2°2((~^)^) I^")- We recall that (j) solves the eikonal 
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equation (3.21) of the previous subsection. Hence, by making use of the asymptotic expansion of h 
and b written above and by collecting the terms with the same order we obtain the following transport 
equations: 

Dtbo= V " \"^ D^bo + aobo, 

\oc\ = l 

D,b..= Y. '-^^DZb^.^a^b.^^Y. '-^^^D^b,^ ^^^"^0, 

(3.24) |a|=l |a|=l |a|=2 



E9"ai(i, 77) 
^ \ " D^b-k + aob-k + f-k, 

l«l=l 

where f^k S ^c°°([-2"2] s-''(R2''))- From Proposition 2.6 and Theorem 2.7 of the previous section we 
deduce the following statement on transport equations with generalised symbols as coefficients and initial 
data. 

Theorem 3.4. The Cauchy problem 

(3.25) Dts = J2 ^""'f ''^ -Pgg + «o(i, X, r])s + f 

\a\ = l 

(3.26) s(0,-,-) = so, 

where 

(i) m€R, f € Qc'='{[-T,T],S"'{K2")) and sq G ^c°°([-T,T],S"'(M2n)), 
(a) ai e ^c~([-T,T],si(M2n)) and ao e ^c°°([-T,T],so(R2n)), 
(Hi) ao is of log-type, 

has a solution s G Sc~([-t,t],5™(k2"))- If we replace Q with Q^'^ in {i), (ii), (Hi) and the assumption of 
log-type with slow scale log-type then s G ^c~([-tt] S"'(m2"))- 

Proof. By applying Theorem 2.7 to the Cauchy problem (3.25) we have the following solution formula: 

s{t, X, rj) = b{t, X, rj) ( so{t, j{t, x,r]),r]) -\- i [ {^'^' ^fi*j^^V,r),r]) 
where, for i = 1, n, 

7i (i , a;, r/) = + / dj^^axia^r]) da, 
Jo 

'yi{t,x,T],T) = Xi + / djj^ai{a,ri) da - / dn^ai{a,vi) da, 
Jo Jo 

b{t,x,r]) = eiIoMsMt,x,n,s),r^)ds_ 

One easily checks that h{t, x, rj) and so{t, 71 (t, x, 77), 72(^5 x, rj),ri) are elements of the generalised symbol 
spaces ^c~([-T,T],so(R2")) and Qc^([-t,t\,S"^{'s?^)), respectively, and that 

f{r,^i^x,^,T)rU) ^ I' J. , 
Jo KT,l[t,x,ri,T),T),r]) J„ 

belongs to Gc°^([-t,t],S"^{V^"))- From this it follows that s(t,x,r]) is a generalised symbol of order m 
satisfying the initial condition s(0, x, rf) = sq{x, rj). □ 
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We can now go back to the system of equations (3.24). With the help of Theorem 3.4 we can solve the 
equations in (3.24): 

Proposition 3.5. If ai and ao are slow scale regular generalised symbols of order 1 and 0, respectively, 

satisfying the assum,ptions of Theorem, 3.1{iii) then there exist symbols b^j G ^c'^([_tt] s-j(R2"))' ^ 
which solve the transport equations (3.24) and such that bo(0,x,r]) = 1 and b-j{0,x,r]) = for all j > 0. 

3.3 Construction of a generalised FIO parametrix and solution formula for 
the Cauchy problem 

We now have all the tools for constructing a generalised FIO parametrix for the Cauchy problem (3.20). 
Combining Theorem 3.1 with Proposition 3.2 and Proposition 3.5 we obtain the following statement. 

Theorem 3.6. Let uq g tJc(IR") and let ai and ao fulfill the assumptions of Theorem 3.1{iii) on ai and 
uq. Then there exists a generalised Fourier integral operator F^{b) with slow scale phase function and 
symbol b of order such that 

DtF^{b)uo = ai{t, Dx)F^{b)uo + ao{t,x, Dx)F^{b)uo + r{t,x, Dx)uo, 
F4,{b)uo{0,-) =uo + ro{t,D^)uo, 

where r G ^^^^|_T^T]^5'-oo^]^2n^^ and ro G Q^-oo^-^2ny 

CoroUciry 3.7. Let the hypotheses of Theorem 3.1{iii) be satisfied. Then the solution u G Q2,2{{—T, T) x 
M") of the Cauchy problem 

Dtu = ai{t, Dx)u + ao{t, x, Dx)u 

u(o,-) = uo eac(M"), 

is equal to F^{b)uo modulo g^2ii-T,T) x K")). 

Proof. From Theorem 3.6 we have that v — F^{b)uo — u£ Q2.2{{—T, T) x R") solves the Cauchy problem 

Dtv = ai{t,Dx)v + ao{t,x,Dx)v + r{t,x,Dx)uo, 
v{0,-) =ro{x,Dx)uo, 

where r{t,x,D:^)uo G G^2ii-T,T) x R") and ro{x,D^)uo G 0^2(I^")- Theorem 3.1 yields the 0°°- 
regularity of v, i.e., v G gk:'2{{-T, T) x R"). ' □ 

4 Microlocal investigation of the solution of a generalised hy- 
perbolic Cauchy problem 

This section provides a microlocal investigation of the solution u G ^(]R"+^) of the hyperbolic Cauchy 

problem studied in Section 2 and Section 3. First wc will concentrate on the microlocal properties of u, 
viewed as a generalised function in both the variables t and x (WFgoou), and second we will fix t and 
investigate the generalised function u{t,-) G ^(M") microlocally {WF goo u{t,-)). Since u can be written 
as the action of a generalised Fourier integral operator F^{b) on the initial data uq G ^c(]R") we focus 
our attention on the microlocal properties of generalised Fourier integral operators. We begin with some 
abstract theoretical results that we will finally apply to the special case of u = F^{b)uo under suitable 
assumptions on the phase function (f). 
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4.1 Microlocal properties of generalised Fourier integral operators: the wave 
front set 



We begin by recalling some results obtained in [13]. Let f2 be an open subset of M". As a preliminary 
step we observe that when (p e x Rp), U CU Q, T CW" \0, V C n xW \0, then 



VGU,ieT + 



■ r \^-^yMy,0)\ 

yeu,ier |^| + \e\ 



is a well-defined element of C. 

In the sequel denotes the set of all points (xo , ^o) & T* (51) \ with the property that for all relatively 
compact open neighborhoods U{xo) of xq, for all open conic neighborhoods r(^o) C R" \ of ^o, for all 
open conic neighborhoods V of conesuppa fl such that V fl {U{xo) x \ 0) 7^ the generalised 
number 



(4.27) Inf 



\^-Vy^iy,0)\ 



yeu{xo),ieT(io) 1^1 + \0\ 

(y,6l)eVn((7(xo)xRf\0) 

is not slow scale-invertible. 

Theorem 4.1. Let I,j,{a) be the functional in C{Qc{fl),C) given by 



[ e'^^y'^'>a{y,i)v{y)dydi, 



where (j) is a slow scale generalised phase function and a is a regular symbol of order m. The Q°°-wave 
front set of I^{a) is contained in the set W|^^. 

The previous theorem can clearly be stated for a slow scale generalised phase function 4>{x, y, ^) (in the 
variable (y,0) ^"^^ ^ regular amplitude a{x,y,^) onQ' xilx MP, with an open subset Q' of IR" . In this 

case 

I^{a)v= [ e''^^'''y'i^a{x,y,Ov{x,y)dxdyd^, 
is the kernel of the generalised Fourier integral operator 

A : g,{n) ^ g{n') -.u^ [ e^^^^'^'^^aix, y, ^)u{y) dy S^. 

We will therefore use the notation Ka for the functional I^{a). The next theorem relates WF go^Au with 
WFq^Ka and WFgoou. 

We recall that given E Cil' x^xR"' x R", E' C (fi' x R"') x (17 x R") is the set of all (x, y, rf) such 
that (.X, y, ^, -ry) eE. A subset E' defines a relation o G of x R" with D.' x R"' when G C 17 x R". 
In detail, 

E'oG:= {{x, e O' X R"' : 3(y, r,) & G, {x, y, r,) G E']. 

Theorem 4.2. Let ^{x, y, ^) be a slow scale generalised phase function in the variable {y, ^) and a{x, y, ^) 
be a regular amplitude on f2' x 17 x R^. Let A be the corresponding Fourier integral operator with kernel 
Ka £ >C(C/c(17' x 17), C) and u a generalised function in Qc{^). If 

(4.28) WF^c» Ka o WFgcx, u C T* (17') \ 
then 

(4.29) WFgoo C {y^F'gooKA o WFgoou) U {WF'gooKA o (singsuppgoo u x {0})), 
where the right-hand side is a conic subset ofT*{fl') \ 0. 
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Proof. Let i be a neighborhood of singsuppgoo u. Let A be the right-hand side of the assertion (4.29) 
with singsuppgoo u replaced by L. It is clear that A is conic and the inclusion A C T*(0') \ follows 
from (4.28). 

Wc will prove that if (.to,^o) ^ A then (.xo.^)) ^ WFgc^ Au. From this it will follow WFg^ Au C A. At 
this point since L is chosen arbitrarily and u has compact siipport wc get the desired inchision (4.29). 

Let Ki be a neighborhood of xq, and K2 = tto (sing supp goo Ka H {Ki x L)). From the assumptions 
on phase function and amplitude we have that A maps Q°° into Q°°. It follows that for a cut-off 
function -0 which is identically 1 on a neighborhood of singsupp^oo u wc can write singsupp^oo Au = 
singsupp^oo A{ipu). If K2 = then taking ip such that supp('(/'u) C L we have that Auj/^j^ can be written 
as / KA{x,y)'ijju{y) dy, where KA{x,y) is ^°°-regular. Hence Au\k-i ^ and (a;o,^) ^ WFgoo^u for 
each ^ ^ 0. 

We assume that K2 is not empty and for j/o € K2 we set 

So = {(C,^): (xo,j/o,e,r/) e WFgooi^^}, 
To = {?? : (yo,??) e WFgoou}. 

We claim that there exist: 

- a conic open neighborhood Eg C R" +" \ of Sq, 

- a conic open neighborhood Fq C R" \ of Fq, 

- a conic neighborhood VF C R" \ of ^0, 

- a number (5 > and neighborhoods Vi and V2 of xq and yo, respectively, 
such that 

(4.30) n (S^ o f 0) = 0, 

(4.31) Fi X ^2 X {(^,?7) : ^ G |?7l < n WFgooi^^ = 0, 

(4.32) So n {(C, 7?) : -7? G To, |^| < 5\r,\} = 0. 

But (a;o,^o) ^ A implies (a;o,Co) ^ WF^oo^l'^ o WFgoou. Therefore, since Sq, Fq and Sq o Fq are 
closed and conic there exist neighborhoods W, Sq and Fq as above satisfying (4.30). If (xo,^o) ^ A then 
(a;o,2/o,€o,0) ^ WFgooJs:^. Indeed, (xo,t/o,Co,0) G WFgooJsT^ implies {xo,^o) G WFgooii:;io(Lx {0}) C A. 
It follows that there exist neighborhoods Vi and V2 as above and 5 > such that (4.31) holds. Finally, 

we have {(0,77) : -7] e Fq} n Sq = 0. Indeed, if (xq, 0, yo, ??) G WFg;ooif^ and {yo,—il) S WFgoow then 
(xo,0) e WF'cjooi^A o WFcjooM. This is absurd since WF^ooi^T^ o WFgoou C T*(0') \ 0. This assertion 
yields (4.32) for a suitable choice of Sq and Fq. By shrinking Vi and V2 wc can also impose 

(4.33) ViCKi, {V2 X (ro)'^) n WFgoou = 0, {Vi XV2X (SoT) n WFgo.KA = 0, 
where the complements are intended in R" \ and R" +" \ 0, respectively. 

Letting a G C^(Vi) identically 1 near xq we have to prove that aAu is ^^Q-regular on a certain conic 
neighborhood of ^o- Let ip G C^(fi) identically 1 on a neighborhood of singsuppgoo u and with supp C 

L. Since aA{{l — tp)u) G G^{fl') we have aAw = aA{ipu) modulo . We write 

a]4Vm(0 = (iV'A,Q!(a;)e-'^^V(y)^i(y))- 
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From K2 = 7rn(sing suppgoo Ka H {Ki x L)) we have that we can restrict our attention to y S K2- Taking 
a partition of unity in a neighborhood of K2 (based on the neighborhoods V2 with the property above) 
we observe that aA'tl)u{^) is a finite sum of terms of the form 

where / G C^iV^ x V2) and j3 G C^(V2)- Applying the Fourier and inverse Fourier transform we can 
write 

{fKA, a{x)e-'^^l3{y)u{y)) = [ fK2{0, r,)a{^ - e)Tu{-r]) SO Sr^. 

JM»'+'' 

Now, from (4.33) we have that the tempered generalised function fKA is Q'^Q-reg^Aax outside Eq and 

that — 9)/3u{—ri) is C/I^Q-rcgular with respect to (9,7]) outside {(6*, 77) : —77 £ Tq, \0\ < S\ri\}. This 
fact combined with (4.32) guarantees that the integral above is absolutely convergent. Note that 



m-e)\<CN{i + \c-e\) 



-N 



for arbitrary N G N and that if W is a sufficiently small conic open neighborhood of ^0 with W C W 
then there exists c> such that \S, - 0\ > c(|^| + \9\) for all £, ^ W and 9 e W. Hence, the estimate 

|S(c-0)|<4(i + |^| + |0|)-^ 

holds for all ^€W' and9 gW. It follows that 

(4.34) / fK2{e, v)a{^ - 0)^{-v) ^0 dr, € Q^oiW). 

Je^w 

In order to understand the integration over W x we begin by proving that fKAi9, r])pu{—r]) is Q'^q- 

regular in {9, rj) for ^^ £ and 77 £ R". From (4.31) we have that JKa is ^^Q-regular in [9, rf) ioi 9 & W 

and |?7| < 5\9\. If jryj > 5\9\ and —rj ^ Fq (4.33) implies that j3u is Q^f^-ieg\\\ai in 77 and therefore also in 

(6', 77). Finally, when 9 and -77 £ Fq then (^,77) ^ Sq from (4.30). Since /Ka is t/^Q-regular outside 

So we conclude that fKA{9, ri)f3u{—r]) has the same regularity property. At this point it is clear that the 
generalised function given by 



(4.35) / fKA{9, ri)a{^ - 0)/3«(-7?) S9 ctr] 

Jeew 

is C7;^Q-rcguIar on W. A combination of (4.34) with (4.35) yields aAu £ Q^fiiW). This means (a;o,^o) ^ 
WF goo Au and completes the proof. □ 

We recall that the ^°°-microsupport of a generalised symbol a on O x M" is the complement of the set of 
points (a;o,^o) with the following property: there exist a representative {ae)^ of a, a relatively compact 
open neighborhood C/ of xq, a conic neighborhood F C M" \ of ^0 and a natural number N G N such 
that 

Vm £ RVa,^ £ N"3c> OB?? £ (0,l]V(a;,0 £ x FVe £ (0,r?] \d^d^a^{x,^)\ < c(C)'"e-^. 
It is denoted by /zsuppgoo a. 

Remark 4.3. Condition (4.28) is satisfied when A = a{x, D) is a generalised pseudodiff'erential operator 
with regular symbol. In this case we have from Remark 5.15 in [13] that if {x,y,^,r]) £ WFgooKA = 
WFgooiir„(j. £)) then x = y and r] = with ^ 7^ and {x,^) £ /xsuppgoo a. Hence 

WF'gooi^a(^,i5)oWFeocuC {(a;,0 : {x,x,^,-0 eWFg^K^^^^n)} nWFg<>.u 

which implies WFgoo/ra(x,D) o WFgoou C T*(0) \ 0. Since 

WF goo Ka{x,D) o (singsuppgoo u x {0}) = 
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Theorem 4.2 yields, for u G Gc{^), 

WF g.oa{x,D)u C WFgooifa(x,D) O WFgcx^M 

C {{x,^) : {x,x,^,-£,) e WF(jo.Xa(x,Z5)} n WFe-u C /xsupp^oo an WFgoow. 
Note that the same inclusion was already obtained in [12, Theorem 3.6]. 

Wc now consider a special class of generalised Fourier integral operators which fulfill (4.28) and for them 
we rephrase Theorem 4.2. This requires some preliminary results concerning the set C<^ of singular points 
and the set W|^„ defined before Theorem 4.1. 

Proposition 4.4. Let [{(j))^] S $(1^xMp) and (p € C^{nxRP\0) such that lim^^o <l)e = <l) in C^{nxRP\0). 
Then 

(t) C G T*(0) \ : V50(x,O = 0}, 

(ii) for any generalised symbol a, 

^[Wela C {(x, V,0(x,0)) -.9^0, {x,9) e conesuppa, Ve(P{x,9) = 0}. 

Proof, (i) From the limit property we obtain that sup,j,g|f \^\=i ~ 01 tends to for any 

compact subset K of fi. If for some {xq, ^o), Vj(/>(xo, ^o) then | V{^(a;, > c on a relatively compact 
neighborhood U{xo) and on a conic neighborhood r(^o)- Since for all n G N there exists £„ G (0, 1] such 
that 

sup |Ve<^e(a;,O-V£0(a;,OI < 
ccec/(xo),?er(€o) 

for all e G (0,£„), we get 

1 



|V5(/.e(a:,OI > |Ve</)(x,C)| - \yiMx,0 - V^0(.t,OI > c- 



n 



for all X G U{xo), ^ G r(^o)) £ € (0,£„). Choosing n large enough we have that the net 
is slow scale invertible, that is {xo,^o) ^ ^[{<i>s)s]' 

(ii) If {xo,^o) S ^[^(0)^] a ^^^^ we find sequences of neighborhoods Un{xo), r„(^o)j and V^(conesuppa fl 
and sequences of points a;„ G C/„(a;o), ^ (^o), G K such that 

|en-V,,(/)eJa;„,^„)| < |loge„r^(|^„| + |^„|), 
where e„ tends to 0. It is not restrictive to assume that 6n has norm 1. Hence, 

(4.36) \^n-^x<l>eAXn,On)\ < | log e^l'^ (1^^ | + 1). 

Passing to subsequences we have that a;„ converges to xq and 9n to a certain with \9o\ = 1 and 
{xo,6o) G conesuppanCj*!?^ -) ]. The first assertion of this proposition implies that V0(p{xo,6o) = 0. The 
sequence |^„| is bounded. Indeed, 

l^nl < ICn - V^^e„(a;„,6l„)| + |V^0e„(a;„,6'„)| < | log£„|-i(|C„| + 1) + \'^x(peAXn,9n)\, 

where, from lim^^tx. |Vx<?i>£„ (xn, ^n)| = \^xH^o,9o)\ it follows that \^x(l>en{^n,Sn)\ < 1 for n large 
enough. Thus, 

|C„|(1- |loge„r^) < |log£„|-i + l 

results in an upper bound for |^„|. As a consequence, passing again to subsequences, we find some ^' = A^O; 
A > 0, such that ^„ — >■ Passing to the limit in (4.36) we deduce ^' = A^o = Va;0(xo, ^?o)- In conclusion, 
we have proved that if {xo,£,o) G ^^[^(0)^] a ^^^^ there exists 6'q ^ such that {xo,9'q) G conesuppa, 
V0^{xo,9'Q)=Oand^o = '^x'P{xo,O'Q). " □ 
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We are ready to consider a Fourier integral operator with phase function 
Prom the previous statements we obtain the following microlocal result. 

Theorem 4.5. Let [{(f>s{x,^) — y^)e] be a slow scale generalised phase function in the variable (y,^) and 
a{x,y,£,) a regular generalised amplitude on fl' x fl x R". Let A be the corresponding Fourier integral 
operator with kernel Ka € C{Qc{^' x r2),C) and u a generalised function in Gc{^)- If there exists 
(j) e C^{Q.' X M" \ 0) with Vx(l){x, ^)y^0 for^j^O such that (j)^ ^ (j) in C^{Q.' x M" \ 0) then 

(i) V^¥'g^KA o WFcjoou C T*{n') \ 0, 

(ii) WFgooAu C {{x,SI^(t>{x,e)) : {x,SI e(t){x,e),e) e conesuppa, (Ve<^(a;, 6*), ^) 6 WFgoou}. 

Proof. An application of Proposition 4.4 to the phase function [(</>£ (a;, ^) — y^)s\ yields 

(4.37) WFgooKA^{{x,y,V:,(i){x,e),-e) : 6^0, {x, y, 6) e conesuppa, Ve<i>{x,e) = y}, 

i.e., 

WFgooKA<^{{x,'^04>{x,e),'^a:4>{x,e),-6) \ O^Q, (x, y, 61) € cone supp a}. 

From the hypothesis on V x(i> wc have that if (.t, ^) G WFgooi^A ° WFgoou then ^ = Vx4>{x, 0) for some 
0^0. It follows that ^ 7^ 0. From (4.37) and Theorem 4.2 we easily obtain the inclusion (ii). □ 

We now apply Theorem 4.5 to the case of a generalised Fourier integral operator solving a hyperbolic 
Cauchy problem as in Section 2. We obtain the following microlocal result. 

Proposition 4.6. Let 

u{t,x)= [ e*l('^-(*'^'''»-l6(t,a;,7?)u5(j?)5-r?, 

where b is a regular generalised symbol of order and [{(j)^ (t, x, rf) — yrj);,] is a slow scale generalised phase 
function in the variable {y, rj). Assume that there exists </> £ C^(IR. x x \ 0) with V(t_j:)0(i, x, r/) ^ 
forr]y^O such that (t)^ ^ cf) in C^K x R" x R" \ 0). Then 

WFgoo?i C {{t,x,'Vt4'{t,x,rj),\/x(t>{t,x,ri)) : {t,x,T]) Gconcsupp6, {'V,^(l){t,x,ri),ri) G WF^ooUo}. 

Example 4.7. We demonstrate how Proposition 4.6 can be used in a simple example. It concerns a 
transport equation with discontinuous coefficients and distributional data taken from in [28, Section 5]. 
In detail, let 

Dtu{t,x) = -H{t-l)D^u{t,x), 
u{0, x) = S{x), 

with t G [0, +oo) and x G R. Here H denotes the Heaviside function and 6 the Dirac measure. Take 
p e 'D(R) with suppp C [—1, 1], nonnegative, symmetric and with integral equal 1. Let {uj~^{e))s be a 
slow scale net with lim^^Q uj{e) = 0. The solution u G G{[0,+oo) x R) of the Cauchy problem can be 
written, at the level of representatives, as 

u,it,x)= [ e'(^-'^^('» Ve(r/)crr?, 

where 

{0, 0<t<l-uj{e), 
Il-u.(e)^eiz)dz l-Cj{6)<t<l + ^{e), 
t-1 t>io{s) + l, 
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and 

r{t-l)/ui{e) 

\S)= / P{z)dz. 

J —CO 

In particular Ae{t) converges to 

'o, 0<t<l, 



m := 



t-l, t>l, 



for £ — > 0. When t ^ 1 we are under the hypotheses of Proposition 4.6. Indeed, y(t,x){x — A{t))r] = 
{-A'{t)r], ?7) 7^ for 77 7^ 0. It follows that 

WFgoo (u|o<t<i) C {(f, 0, 0, ?7) : ?7 ^ 0, < t < 1} 

and 

WF5c» {u\t>i) C {{t, t-l,-v,v)- Vt^O, t> 1}. 

The fact that A is not difFcrcntiablc at t = 1 docs not allow us to apply Proposition 4.6 in a neighborhood 
oi t = 1. The wavefront set at t = 1 has been computed by direct calculations in [28]. 

Another class of generalised Fourier integral operators satisfying condition (4.28) can be found in solving 
a hyperbolic Cauchy problem where the principal part depends only on t. More precisely, we will deal 
with phase functions of the type 

[{x'n + ipe{t,r]) -yr])s], 

where ip^ is homogeneous of order 1 in 77. Note that differently from Proposition 4.6 we do not require 
any convergence as £ — > but simply that ipe does not depend on x. Let us consider the Cauchy problem 

n 

Dtu = ^ aij{t)DjU + ao{t, x)u, u{0, •) = uq e Gc^W)- 

i=i 

A typical example of <^ = [((/Je)^] is given by 

" / /■* N 



'7) = ^ oij (a) do- 



where 'yj{x,t,s) = Xj + a-ij{a)da are the components of the characteristic curves corresponding to 
the operator Dt — ^2^=1 

Proposition 4.8. Let 

u{t, x)^ [ e*[(^''+'^= 6(t, X, v)Mv) ■= Buq, 

where b is a regular generalised symbol of order and [{xrj + ip^{t,ri) — yr])g] is a slow scale generalised 
phase function in the variable {y, rj). If the net (sup|j|<y ^_^q dt^fsit, ^))£ ^■s bounded for all T > then 

(4.38) WFgoou C {WF goo Kb o WFgooWo) U {WFgoo Kb o (singsuppgoc uq x {0})). 

Proof. The microlocal result (4.38) is directly obtained from Theorem 4.2 by checking that the inclusion 

WF^ooKboWFq^uo C T*(K X R") \ 

holds. We will prove something stronger. We will prove that if {to,xo,yo,To,^o,ilo) G WFgoo then 
Co = -m- Thus (to,a;o,yo,0,0,?7o) ^ WFgooKB if rjo 0. 

By Theorem 4.1, for the point (toi 2^0, yo; "n], Co: Jyo) G WFgoo ii'^ we can find sequences of neighborhoods 
U„{to-xo,yo), r„(ro,Co,%) and sequences of points {tn,Xn,yn) & Un{to,xo,yo), {Tn,^n,Vn) & r„ and 6'„ 
such that 

)|<|log£„| ^(|(-r„,Cn,77„)| + |^„|), 
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where e„ tends to 0. Again we may assiimc that 0„ has norm 1. Hence, 

\{Tn,^n,Vn) " (9t<(Ce„ (f„, 6'„), 6'„, -6'„)| < | loge„|"^(|(T„, + 1). 

Prom the hypothesis on dt^Pe it follows that 

\{rn,^n,Vn)\ < l^OgSnl''^ {\{Tn, ^n,Vn)\ + 'i-) + \dtV>Sn{tn,6n), On, -On)\ < | logEnl"^ |(t„, ,^„, 77„)| + C, 

which means that the sequence | (t„ , ^„ , ?7„ ) | is bounded. Passing to subsequences we have that {Tn,£,n,Vn) 
converges to A(to,^0)?7o) and On to some 0. Finally, 

\{^n,Vn) - {On,-On)\ < \{Tn, ^n,Vn) - {dtfcnitn, On), On, -On)\ < l^OgEn]''^ {\{Tn, ^n,Vn)\ + 'i-) 

yields 

{X^o,Xvo) = iOo,-Oo), 

i.e., ,^0 = □ 
4.2 The Hamiltonian flow and the ^°°-wave front set at fixed time 

We now fix t and investigate the microfocal properties of the solution u{t, •) of the Cauchy problem 
(2.10). As in the previous subsection we begin with some theoretical results on generalised Fourier 
integral operators of the type 



and we then consider the specific case of a strictly hyperbolic Cauchy problem. The following theorem is 
modelled on [23, Theorem 4.1.6] and makes use of the composition formula for generalised operators in 
[10]. 

Theorem 4.9. Let f2 be an open subset of R", (j) a slow scale regular generalised symbol on x RJJ 
homogeneous of order \ in rj, b a slow scale regular generalised symbol of order on fl x R" compactly 
supported in x and u € Gc{^)- Let (xq, ^o) G T*{^) \ and {yo,ilo) ^ WFgoou. // there exist 

- a representative {4>e)e of (j) such that Vx^e and 

( inf \s/,cj>,{x,^)\)^ 

is slow scale strictly non-zero for all K (sV,, 

- conic neighborhoods W and Z o/(j/o,??o) o,nd {xo,^o), respectively, with W C (WFgoouy, 

- £o G (0, 1] and a slow scale net {Xs)s 

such that the following implication holds 

(4.39) {x,V:cMx,v)) & Z V^(/>e(x,77)| > AJ^ for all e G {0,eo] and all [y,?]) ^W, 
then 

(4.40) {xo,^o) ^ WFgoc {F^{b)u). 

Proof. By definition of the C/°°-wavefront set we can find a conic neighborhood Wi of {yo,r]o) such that 
W C Wi C (WFgoou)'^. From Lemma 3.4 in [12] we can find a symbol c e S^{n x R") with /xsupp cCWi 
and such that c{x, ri) = 1 on W when \r]\ > 1. In addition, we can construct c so that supp^. c (s fi. From 



23 



Theorem 3.6 in [12] we have WFgoac{x,D)u C WFgoou n /isupp c = 0. It follows that c{x,D)u G g^{n). 
We now write 

F4b)u = F^{h)c{y, D)u + F^{b)e{y, D)u, 

where e{y, rj) — l—c{y, rj). Since c{y, D)u E Q^{i}) we obtain from the mapping properties of a generalised 
Fourier integral operator that F^{b)c{y, D)u is t/°°-regular as well. This means 

WFgo.F^(6)u = WFeooF^(6)e(2/, D)u. 

We will therefore prove that the implication 

(4.41) (j/o,%)^WFac»(«) ^ (xo,6) ^WF6oo(F^(6)e(y,Z))u) 

is true. Making use of the characterisation of the ^7°°-wave front set stated in [12, Theorem 3.11] it is suffi- 
cient to prove that for some pseudodifferential operator a(x, D) with a S x R") such that supp a C Z 
and a{x, ^) = 1 on some smaller neighborhood Zi when |^| > 1, we have a{x, D)F^{b){e{y, D)u) e Q°°{Q,). 
We are under the hypotheses of the composition theorem in [10] (Theorem 5.11 in [10] without depen- 
dence in t). Hence, modulo a 5°° -regularising operator which does not aifect the wavefront set, we have 
that the operator a{x, D)F^{b){e{y, D)u) is the generalised Fourier integral operator 



F^{bi){{e{y,D)u)){x)= [ e'^^^'^^b,{x,7^) / e'^y^^'eiy, D)u{y) dy drj 



where bi G ^so(qxM") following asymptotic expansion: 

a 

At the level of representatives it means 

a 

and since we argue modulo we can consider bi^e as a sum of a convergent series obtained from the 
asymptotic expansion above as in [10, Theorem 2.2]. Hence bi^^{x,ri) = if {x,'Vx(f>s{x,ri)) ^ Z and 
£ G (0, £o]. In the sequel we will complete the proof by showing that F^{bi)e{y, D) is an integral operator 
with kernel in Q°°{Q, x O). This will follow from the fact that the transposed operator *e(j/, D) *F^{bi) has 
kernel in Q°°{Q x ft). Working at the level of representatives we have that the kernel of *e(j/, D) *F^{bi) 
is given by 

(4.42) f e'^''''^'''''^-y'''>{^e){y,-r])bi,e{x,r])v{x,y)dxdySr], 

JnxQxM." 

where *e{y,i]) ~ X]Q(Q^0~^-O^^y ^(j/, — 77). Again arguing modulo it is not a restriction to assume that 
*e(y, —77) = when {y, rj) G W. It follows that in the integral (4.42) we may assume that (a;, Vx^e{x, r])) G 
Z and {y, rj) ioi e € (0, eoj. Thus, by integrating by parts we have 

[ Qi{M^,v)-yv) ( tg) -rj)bi^^{x, r])v{x, y) dx dy drj 



dr)v{x,y) dxdy, 



where we can make use of the hypothesis (4.39). Since the nets involved in the integral 

(4.43) / e^(..(.,.)-..) ^^0-<^;-^)^^^-(^;;)) ^, 

are of slow scale type and |Va;(/'e(a;, rj) — y\'^^ > X^^^ where {Xe)e is a slow scale net as well, we conclude 
that (4.43) generates a generalised function in g°°{fl x O). □ 



24 



Theorem 4.9 applies to the sohition 

u{t,x) = F^{b){uo){t,x)= I e'^^*'^'^^b{t,x)^{r,)dr] 



(sec Proposition 2.6) of the Cauchy problem (2.10) under a certain convergence assumption on the 

Hamiltonian flow. 

Theorem 4.10. Let u G 5(M"+^) be the unique solution of the Cauchy problem 

n 

DfU = ^ 0-1, jit-: x)DjU + ao{t,x)u, 

u(0,-) = wo eac(R")- 

Under the hypotheses (/i2) on the coefficients aij, j = l,...,n, and {i2) on the coefficients ao, let 

(j){t,x,ri) — '^"i^i'yh{3:,t,0)rih be the solution of the eikonal equation (2.11) and b{t,x) G Q°°(W^~^^) 
as in (2.15). Let 4> have a representative {4'e)e satisfying the following two conditions for any t G M.' 

(i) Vx<^e(i, x,ri) ^ for all e, x and ri ^0 and the net 

inf \Va^<j)e(t,x,-P-)\ 

is slow scale strictly non-zero; 

(ii) the Jacobian {dj'yi^e{x,t,0))ij is invertible and the solution {y,ri) = Xt,£{x-,^) of the system 

£, = yx<t>e{t,x,r]), 
y = \7^(f)^{t,x,r]), 
converges to a limit Xt in C(T*(M") \ 0,T*(IR") \ 0) as s ^ 0; 
(Hi) Xt defines a bijection on T*(M") \0. 

Then 

WFgo.u{t, •) C Xt"^(WF5ocUo). 

Proof. We denote the linear map corresponding to the Jacobian {dj'yi^s{x,t,0))ij by rs{t,x). The in- 
vertibility of r^{t,x) allows us to solve the Hamilton- Jacobi system above with respect to the variable 
{y,r]). More precisely we have 

Xt,e{x,0 = {H,6{x,t,0), ...,-yn,eix,t,0), (F* )"^ (i, a;)^ 

and 

Xm (2/' ^) = ("Ti.^ (y> 0, t), 7n,e(y, 0, t); r* {t, 7^ (y, 0, t))r]), 
where 7e(y,0,i) = (7i,e(y,0,i),...,7„,e(y,0,i)). 

Let (.To,^o) G T*(IR") \ and (j/o,f/o) = Xi{xo,£.q)- The limit in (ii) entails the following assertion: for 
all conic neighborhoods W of (j/0)'7o) there exists a conic neighborhood Z of (a;o,^o) and some Sq G (0, 1] 
such that Xt,e{x,0 G W for all (x,^) G Z and e G (0,eo]- In other words, if {x,Wx(l>£{t,x,r])) G Z and 
e G (0,eo], then (V^(/)e(f, a;, ry), r/) G W. Thus, {x,Wx<Ps{t,x,r])) G Z and s G (0,£o] imply for all Wi with 
W C Wi the assertion 

3c>0y{y,r])^Wi \y -yr,Mt,x,v)\ > c 

Combining the hypothesis (i) with this result we arc under the assumptions of Theorem 4.9 for fixed t. 
Therefore, (t/o,%) = Xtixo,^o) ^ WF^oo-uo implies (.to,^o) ^ WFg;ooM(t, •) or in other words 

(4.44) WFgo.{u{t, •)) C xT'O^Fg^Uo). 

□ 

Remark 4.11. When the coefficients aij and ao are C°°-functions Theorem 4.10 recovers the propagation 
of singularities result obtained in [12, Proposition 4.3]. In this case arguing by time reversal one proves 
the inclusion x7^(WFgooUo) C WFgoo«(t, •) and therefore WPgoou^t, •) = Xt~^(WFgooUo)- 
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4.3 Examples 

This subsection contains some examples of first order hyperbolic Cauchy problems whose principal part 
depends only on t. Applying the results of the previous sections we can express the Colombeau solution 
by means of a generalised FIO formula and perform a microlocal investigation. 

Proposition 4.12. Let u G ^(M"+^) be the unique solution of the Cauchy problem 

n 

DfU = ^ aij{t)DjU + ao{t,x)u, 

«(o,-) = «o eac(]R"). 

Under the hypotheses (/i2) and (z2) let 4>{t,x,r]) = Y^'l^iJh{x,t,0)r]h be the solution of the eikonal 
equation (2.11) and let b{t,x) G Q°°{M."'+^) be as in (2.15). If there exists a choice of representatives 
{(ii,j,s)e such that Aj^s{t) = — /J aij^e{z) dz, j = 1, ...,n depends continuously on e and 

lim Aj^s{t) = Aj(i), for all i G M, j = 1, n, 

then 

WFaoou(i,.) = Xt"'(WFaoo«o), 

where 

XtHv, V) = {yi + Mt), Vn + An{t),Tl). 
Proof. The components 'yj,s{x, t, s) of the characteristic curves are the solutions of the Cauchy problems 

^7j,e(a;,i,s) = -aij,£(s), ^j^e{x,t,t) = Xj. 

Hence, 

jj^s{x,t,s) = Aj^s{s) + Xj - Aj^^{t) 

and 

n 

(t>e{t,x,ri) = ^{xj - Aj^^{t))r]j. 
It follows that the hypotheses (i) and (ii) of Theorem 4.10 are satisfied with Xt,e{x,^) given by 

y = le{x,t,0) = X-As{t). 

Now 

Xt{x,0 = limxt,c(a;,0 = (x - A{t),0 

is a bijection on T*(K") \ with Xt^iv^v) = {v + A{t),rj). From Theorem 4.10 we have 

WFg^{u{t,-)) C Xt"'(WFgooUo) = {{y + A{t),r]) : {y,r]) G WFg^uo}. 

We will now argue by time reversal. We fix to S M and we set v{t,x) = u{to — t,x). From the original 
Cauchy problem we obtain 

n 

DfV = — ^^oi,j(io — t)Djv{t, x) — ao{to — t, x)v 
with initial condition v{0,x) — u{to,x). The corresponding characteristic curves are 

^j,e{x, t, s) = Xj+ / aj^e{to - z)dz- / aj,e(io - z) dz 
Jo Jo 
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and therefore 

7j,e(a;, t, 0) = Xj - Aj,e(to -t) + AjAto)- 
In particular, denoting the corresponding Hamiltonian flow by Xto^e^ have 

xtoA^^O = {x + ^e{to),0, Xto]e{y^v) = {y - ^eito),v)- 

Moreover 

Xto{x,0 ■= limXto,e(a;,0 = lim(a; + A£(to),0 = {x + Hto),0 

e— >0 £— >0 

and 
Hence, 



or in other words 
which imphes 
Since Xto = Xto deduce 
In conclusion. 



WFe=»(t;(io, •)) C XtV(WFe=»(^;(0, •))) 

WFe<»K)cxtV(WFe»(«(to,-))) 
Xto(WF5«(uo)) C WF5oc(u(fo,-))- 
X4V(WF5»«o) C WF6oo(u(to, •))• 
WFg^w(t,.) = xr'(WFe^«o). 



□ 



Example 4.13. A transport equation with discontinuous coefficient and distributional data. 

We continue studying the Cauchy problem from Example 4.7. By applying the previous proposition we 
can undertake a microfocal investigation of the solution u £ Q{{0, +oo) x R). Using the nets and notations 
of Example 4.7 we easily see that the assumptions of Proposition 4.12 are fulfllled. In particular, the 
Hamiltonian flow is given by 

xt,e(a;,0 = {x-A,{t),0, xll{y,n) = {y + As{t),v), 

with Ae(t) deflned in Example 4.7 and having limit 



A(i) 

Hence 



0, 0<f<l, 
t-1, t>l, 



Xt{x,0 ■= limXt,£(a;,0 = ix-A{t),^), Xt ^{y,n) = lim Xt,^^ (y, ??) = {y + A{t),r]). 

Concluding, for each t > we have 

WFg^{u{t, •)) = xi'mo-m = {(A(t), r/) : 77 ^ 0}. 
Example 4.14. An example in higher dimension. 

We can adapt the previous example to higher space dimensions. For instance, let 

Dtu{t, xi,X2) = -aiH{t - l)Dx:iu{t, xi,X2) - a2H{t - \)Dx^u(t, xi,X2), 

U{0, X) = tlj{xi,X2)6c{xi,X2), 

where ai, 02 S M, 6c is the Dirac measure along a smooth, simple curve C G M?, and ip G 2?(M^). In this 
case the components of the characteristic curves are given by 

7i,£(a;i,i, s) = aiAs{s) + xi - aiAe(i), 'y2,e{x2, t,s) = a2As{s) + X2 - a2Ae(i) 
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with Ag as above. This yields 

Xt,e{x,0 = {{xi - aiAs{t),X2 - a2AE{t)), (6,6)), Xt,l{y^'n) = {{yi+aiAs{t),y2+a2As{t)), {r]i,m)) 
and the limit 

Xr^(y, V) = {{yi + aiA{t), y2 + a2A{t)), (r/i, 772)). 

In conclusion, 

WFa»Kt,.)) =xr'(WFc;^(#c)) 

= {{xi + aiA{t),X2 + a2A{t),^i,^2) ■■ {xi , X2) e supp ip , (6,6) € iV(xi,a;2)} 

where N{xi,X2) denotes the set of nonzero conormal directions of C. 

Example 4.15. An example with the Dirac measure as a coefficient. 

We consider the Cauchy problem 

Dtu{t,x) = ^S{t~l)D^u{t,x), 
u{Q,x) = uo{x), 

with uq S ^c(K)- In the rcgularisation of S{t — 1) via convolution with a mollifier we use p £ !>(]&) as 
in Example 4.7 and a slow scale net {lo~^{s))s with < oj{e) < 1 and lim^^o Lo{e) = 0. The solution is 
represented by 

Ue{t,x) = Uoe{x - Xe{t)) 

with Xe{t) as in Example 4.7. Since suppp C [—1,1] and J^^p{z)dz = 1/2 we have that Xs{^) = \ for 
all £ € (0, 1] and thus 

'1, t>l. 



lim Ae(rt =: Xit) = < 

s—>0 



1 

2' 

0, t<l. 



2' ^ — 1, 



By Proposition 4.12, 

WFg» {u{t, •)) = {(x + A(t), 6 : {x, 6 e WFg~ wo} 

for all t gR. Note that the limits of the flows Xt,£{x,^) are homeomorphisms at each fixed t, while the 
limiting two-dimensional characteristic coordinate change {t, x) {t,x — X{t)) is still bijective, but no 
longer continuous. Such a situation is admitted by Proposition 4.12. 
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